Abstract 



The equations for various spin particles with oscillator-like interactions are discussed in 
this talk. Contents: 1. Comment on "The Klein-Gordon Oscillator"; 2. The Dirac 
oscillator in quaternion form; 3. The Dirac-Dowker oscillator; 4. The Weinberg oscillator; 
5. Note on the two-body Dirac oscillator. 
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1 Comment on "The Klein-Gordon Oscillator" 

[By S. Bruce and P. Minning, Nuovo Cim. 106A (1993) 711] 

Concept of relativistic harmonic oscillator with intrinsic spin structure has been proposed long ago [Q. 
However, in connection with the publications of Moshinsky et al. [0 the interest in this simple model with 
a j = 1/2 Hamiltonian that is linear in both momenta and coordinates has grown recently ft]. Analogous 
type of interaction has been considered for the case of j = and j = 1 Duffin-Kemmer field 0] and for 
the case of j = Klein-Gordon field §. 

In the paper || the operators Q, coordinate, and P, momentum, have been represented in n ® n 
matrix form 

Q = fjq, P = fjf, (1-1) 
with f/ 2 = 1. The interaction in the Klein-Gordon equation has been introduced in the following way: 



P -> P - im^n ■ Q, 



(1.2) 



where for the sake of completeness fl is chosen by 3 ® 3 matrix with coefficients Ojj = LOiSij (the case of 
anisotropic oscillator). The 7 matrix obeys the following anticommutation relations: 



{7,^} = 0, 
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1. 



(1.3) 



The Klein-Gordon equation for ^(q, t), the wave function which could be expanded in two-component 
form, is then 



o f \ 

— %(q, t) = [p 2 + m 2 q- O 2 • q + mj trVL + m 2 ) $(q, t), 



dt 2 

what gives the energy spectrum p] 

i 2 = 2m (uiNi + <jJ 2 N 2 + u 3 N 3 ) . 



T? 2 



N 1: N 2 , JV 3 = 0,1,2. 



E (b)m- m2 = 2m{cu 1 (N 1 + l)+u; 2 (N 2 + l)+u; 3 (N 3 + l)), 



(1.4) 



(1.5) 



It becomes, up to an additive constant, the spectrum of the anisotropic oscillator in the non-relativistic 
limit. 

However, the physical sense of implementing the matrices r) and 7 in |3j is obscure. In this Section 
we try to attach some physical foundations to this procedure. 

It is well-known some ways to recast the Klein-Gordon equation in the Hamiltonian form [|6|- [ict| . First 
of all, the Klein-Gordon equation could be re-written to the system of two coupled equations [8, p. 98] 



dx a 



dx a 



(1.6) 



where n — mc/h (in the following we use the system where c = % = 1). By means of redefining the 
components they are easy to present in the matrix Hamiltonian form (cf. with |ll|] ) 
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(1.7) 



provided that 



id t ^ + to* 
= -t\/i * =Pi^. 



(1.8) 
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Using matrices 
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(1.9) 



(i, j, k are the orth- vectors of the Euclidean basis) and substituting analogously [2a] , i. e. p ^ p — imu>(3r, 
we come to the equation for upper component 



(E 2 - m 2 



r ->2 , 2 2 ^2 

\p + m to r 



3mw] (/> 



(1.10) 



what coincides with Eq. (10a) of ref. g in the case uj\ = uj 2 = <^3- 

The similar formulation also originated from the Duffin-Kemmer approach [Q . In this case the wave 
function $ = column(<f)i, 4>2, Xi) is five-dimensional and its components are 



0i = (idt^S + m^)/V2 
4> 2 = (id t V - m*)/\/2 



It satisfies the equation 



5,, 



(1.11) 



(1.12) 



(our choice of 5 5 dimension /3-matrices corresponds to ref. [4b]). As shown there, the substitution 
p — ► p — imajr/Qr, where 










= 




V 


— H-3®3 



(1.13) 



leads to the equation ( l.lCj ) for both (fix and fa. Let us remark, in both the approach based on Eq. (1.7) 
and the Duffin-Kemmer approach, Eq. ( 1.12| ), we have the equation 



(E 2 — m 2 )xi = {pi — imujXi)(pj + imuJXj)xj 



(1.14) 



for down component, which seems to not be reduced to oscillator-like equation. 

Then, Sakata-Taketani approach [pi fLOl is characterized by the equation which we write in the form: 



p{r 3 + ir 2 )p 
2m 



771T3 > $ 



(1.15) 



with n being the Pauli matrices. $ = column((j), %) is the two-component wave function with components 
which could be written as following: 



4> = (* + ^*)/x/2 
X = (* - ^9 t *)/V2. 



(1.16) 



From the previous experience we learned that in order to get the oscillator-like equation we need to 
do substitution with matrix which anticommutes with matrix structure of the momentum part of the 
equation. In our case the matrix which has this property is t\ matrix. Therefore, we do the substitution 



p — > p — im<jjT\r 



and come to 



E 2 ^ = [p 2 + m 2 u 2 r 2 - 3muj + m 2 ] £, 



(1.17) 
(1.18) 
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where £ = <f> + \ and to the analogous equation for 77 = <j) — x = — {<t> + x)- 111 the process of calculations 
we convinced ourselves that the interaction Hamiltonian 



H = (p 2 + m 2 oj 2 x 2 — 3muj) (t 3 + ir 2 ) + mr 3 

2m v ' 



(1.19) 



is the same as in [4b, formula (3.9)] since ti(t 3 + 2t 2 )ti = — (r 3 + iTi) and (r 3 + it2)ti = t 3 + iT2. 

Let us also recall the Dirac oscillator in (1 + 1) dimensions [3j]. Curiously, the formula (1.7) in (1 + 1) 
dimensions looks like the Dirac equation in the case of the choice of 7-matrices as in [3j], i. e. 



a = t x (3 = t z 

In the case of the choice of 7-matrices as in Jl|], i. e. 

7o = n 71 = ir 3 , 



(1.20) 
(1.21) 



in order to obtain the Dirac oscillator it is necessary to do substitution (1.17). 



2 The Dirac oscillator in quaternion form 

The quaternion (and conjugated to it) with real coefficient is defined as 

q = Qa + iqi + jq2 + kq 3 

q = qo-m-jq2-kq 3 . (2.1) 

The basis vectors satisfy the equations i 2 = j = k 2 = — 1 and ij = —ji = k with cyclic permutations. 

Considering a two-component^maternionic spinor (or SL(2, H) spinor) one could write the free Dirac 
equation as following, ref. [14c,d]cl 

f • 9* - mr 3 *fc = 0. (2.4) 

Anticommutation relations for T are given in [14d,p.222]. In Pauli representation (i — * — \/— In, j — -> 
— \J — IT2 and k — > — \ / — lr 3 ) it goes through to usual Dirac equation and its complex conjugate. As 
mentioned in |l4f| it is convenient to diagonalize the matrices entering in Eq. (2.4) using matrix 

T-if: (2.5) 



In such a way we come to biquatcrnionic formulation (qi G C, i. e. the coefficients are the complex 
numbers): 

8il> L + im{p f R = (2.6) 
<9t/^ + im^> L = 0, (2.7) 

where ipL = ">jjp+, tpR = ^V-- This decomposition of \& into left ideals is carried out by means of the 
projection operators p± = (bo ± b 3 )/2. New basis is 60 = 1, 61 = \/~^li,&2 = %/— Tj, 63 = -\/— Tfc and 



"Let us mention other proposed equations, e. g. Jl5|]-[|li||. However, the equation given by J. Soufiek [l6[ 

+ jd 2 + kd 3 )ip = (-V^Tr 2 d + r 3 m)V> (2.2) 



For a gauge theory, involving quaternion- valued fields, sec |H|-pl|. 
|l4] the norm N(q) = qq = qq is not well defined. However, m the 



leads to the tachionic solutions of the Dirac equation. 
Meantime, in the biquaternion formulation of Morita 

case when components of quaternionic spinor are the grassmanian numbers, the norm vanishes: 

N{q) = -rfegV = (go) 2 - (91) 2 - (<? 2 ) 2 - (q 3 ) 2 = 0. (2.3) 
See footnote on p. 1651 of ref. [14a] where the necessity of this suggestion is discussed. 
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&o = bo, b a = — b a . Introducing interaction in the form di — * di + TsVi(x), V is the compensating field for 
this type of Sp(l,Q) transformations, and taking into account that the vectors of biquaternionic basis 
anticommute b a bp + bpb a = —2r] a p, r] a p = diag(—l, 1, 1, 1), we come to the equations for the left and 
right spinor-quaternions in the following form: 



(E 2 - m 2 )tp L = [(p 2 + k 2 x 2 ) -3k- 2e ijk b k x i p 3 ] ip L 
(E 2 -m 2 )i> R = [(p 2 + k 2 x 2 ) + 3k + 2e ijk b k x i p j ]ip R 



(2.8) 
(2.9) 



if we choose Vi(x) = kxi. Eqs. (2.8) and (2.9) are the Dirac oscillator equations in the Pauli rep, b k — > T k . 
Analogous equations for ipR and ^r L could be obtained if we start from (2.6) and (2.7) with the opposite 
signs at the mass terms. In the above, we assume, following for Morita ]14|, that the imaginary unit s/—l 
commutes with Hamilton's basis vectors k, what is not obviously. Consideration of y/—l on an equal 
footing with the Hamilton's units would also be interesting. 

The investigation of interactions of quaternionic Dirac field deserves further elaboration. 



3 The Dirac-Dowker oscillator 

In this Section we start from the equation for any spin given by Dirac p3], see also p3| , in the form 
written down by Corson, ref. [8, p. 154], (here we use Corson's notation) 



d AB v B (k + + |, I - |) - m (^±i) 1/2 v A (l)iP(k, 0=0 

d AB v A (l)ip(k, I) + m (^J v B (k + \)^j{k + i, I - i) = 



where va and are the rectangular spinor-matrices of 2k rows and 2k+ 1 columns (see, e. g., section 17b 
of ref. Q for details). The wave function -0(fc, I) belongs to the (fc, I) representation of the homogeneous 
Lorentz group. The choice 1 = 1/2 and k = j — 1/2, j is the spin of a particle, permits one to reduce a 
number of subsidiary conditions. Moreover, the equations (3.1) are shown by Dowker |m| to recast to 
the matrix form which is similar to the well-known Dirac equation for j = 1/2 particle 

a^d^ = ml, 

a^T = (3.2) 

The 4j- component function <i> could be identified with the wave function in (j, 0)©(j — 1, 0) representation. 
Then, T, which also has 4j components, is written down 

T = (-i)^( 2i )-i (^ ( \" ^^/^Vci - i i)- 0.3) 

V u A (j)®v A (±) J 2 2 

and it belongs to (j — 1/2, 1/2) representation. The matrices and = obeys the anticommutation 
relations of Pauli matrices 

a^a u) = g^. (3.4) 

This set of matrices has been investigated in details in ref. [24b, c] and a M was proved there to satisfy all 
the algebraic relations of the Pauli matrices except for completeness. 
Defining p^ = —id^ and the analogs of 7- matrices as following: 

7"=(° ~T) (3-5) 
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the set of equations (3.2) is written down to the form of the Dirac equation 



(p M 7 /1 - m) 



0. 



(3.6) 



However, let us not forget that $ and T are 2-spinors only in the case of j = 1/2. 

As mention ed i n, e. g., ref. [8, p. 33, 124], in the case of spin j = 1/2 the set of 7- matrices in 
representation (3.5) 





ir 2 

IT 2 









\ IT 
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IT 3 


IT 3 






IT 





(3.7) 



is defined up to the unitary transformation and Eq. (3J3) could be recast to the Hamiltonian form given 
by Dirac (with a% and [3 matrices) by means of the unitary matrix. It is easy to carry out the same 
procedure (a k = <S7°7 <S _1 and (3 = S^S^ 1 ) for 7 matrices (Eq. 3^) and functions of arbitrary spin 
(^ = iS _1 <l>). For our aims it is convenient to chose the unitary matrix as following: 



\[2 \i^4j®4j ^4j®4j 

After standard substitution p — > p — imujj r we obtain 

J E<j) — —i [ao(dp) + imoj(df)] v + mao(f>, 
\ Ev = i [ao(ap) — imu>(af)] cf> — mceov. 



(3.8) 



(3.9) 



Since it follows from the anticommutation relations ( J3.4[ ) that ona^ = aocti we have the equations which 
coincide with Eq. (8) of ref. [2a] or Eqs. (3.6) and (3.12) of ref. [2h] except for 77 — > a M , i. e. their 
explicit forms, 



(E 2 



P 



t2 — 3aomuj — mujaoa^a^^sj , 



(E 2 + m 2 )v = p 2 + rncur 2 + 3aomu! + rnu>aoa^a^r l \/ j 



(3.10) 
(3.11) 



Thus, we convinced ourselves that we got the same oscillator-like interaction and the similar spectrum 
as for the case of j = 1/2 particles in [2a]. 



4 The Weinberg oscillator 

The principal equation of 2(2 j + 1)- component approach |25j| in the case of spin j = 1 is 

(l^V^Vu + M 2 )^ =1 \x) = 0, (4.1) 
with 7 a /3 being defined by the formulae 






Sij SiSj 


&ij SiSj SjSi 





( iSA 


744 ^ L 


{ -%Si J ' 





-n='u= { SU, '„•' )• -11^ ( H ® 3 ) (4-3) 



3<g>3 



(Si are the spin matrices for a vector particle). 

The j = 1 Hamiltonian has been given in refs. [^6[ [27]: 



2£ 2 



2^2 _ M 2 



2£ 



2 £2 _ M 2 



(ap)' 



(4.4) 



where 







l3<8>3 





Though this way of description is a little bit antique, attention has again been paid to it recently, 
e. g. |28j| ^9). This formalism presents oneself the example of the Bargmann-Wightman-Wigner type 
quantum field theory [28b]. The remarkable feature is the fact that boson and its antiboson have the 
opposite parities. 

In general, the upper and down components of 6- component wave function do not uncouple neither 
under the interaction p — ► p — imuj(3r nor under ^^^yU^Xu. However, if we introduce the Dirac oscillator 
interaction so that the conditions of the longitudity of ^ = column(<f>i, Xi) respective to f 



rx(p = 0, rxx = 
are fulfilled, we come to more simple equations (£ = <j) — %, rj = <p + x) 



(Sp) + k(Sr) (Sp)?] 



(2E 2 -M 2 )£ = E(Sp) V + (Sp)-k(Sr) (Sp)£, 

E(sm - 

which could be uncoupled to the following form (k = imco) 

(Sp) (E 2 - M 2 )(Sp)£, = (Sp) \p 2 + muj 2 r 2 + Zmuo + AmuoS[r x p\\ (Sp)£ 



(4.5) 

(4.6) 
(4.7) 



(Sp)(E 2 - M 2 )(Sp) V 



(Sp) 



->2 2 -»2 

p + muj r 



Zmoj — Amu)S[r x p\ (Sp)rj 



-imuj(2E 2 - M 2 )(Sr)(Sp)ri 



(4.8) 



(4.9) 



These equations can be considered as the extension of the equations with Dirac oscillator interaction 
to the j = 1 case, for the components (Sp)£ and (Sp)rj. However, remark that one has the additional 
spin-orbit term acting as earlier at rj. 



5 Note on the two-body Dirac oscillator 

The two-body Dirac Hamiltonian with oscillator- like interaction is given by (see, e. g., refs. [2d,f]) 
( 9 d \ 



1-1 i 

— (di + d 2 ) ■ p + ~ "2) • p — ~^2^° ti ~ " 2 ) ' ^ B + m (Pi + #2) 



1>, 



where the matrices are given by the direct products 



Oil 



ax 

(7 1 



1 2 ®2 
1 2 ®2 



a 2 



L2®2 1 

1 2 



<7 2 

02 



(5.1) 



(5.2) 



. _ / II 2 J \ / 1 28 J 



If we are in the center of mass system (c.m.s.) it is possible to equate P = 0. 

Now we apply the same procedure like that was used for transformation the Bargmann-Wigner equa- 
tion to the Proca equations (see, e. g., p. 30-31]). The 16- component wave function of two-body 
Dirac equation could be expanded on the complete set of matrices: (-f^C), (a^C) and C, (7 5 C), and 
(7 5 7 M C). We consider the system multiplied by C, the matrix of charge conjugation, in order to trace for 
symmetric properties under oscillator-like potentials, see also 30, p. 31]. The wave function is decomposed 
in symmetric and antisymmetric parts using the above-mentioned complete system of matrices 

ip = i>{ a f3} +i>[afl, (5-5) 

where 

= ^ v CvP A » + <v C V0 F ^ (5-6) 

-0[a/3] = Ca^ + ^C^ + ^^C^Ap. (5.7) 

In such a way we obtain the set of equations: 

EA = 0, EA = -2m4> (5.8) 

E<\> = 2iV2{fi - if^F 10 (5.9) 

E$ = -2mh + V2e ijk (pi + ir l )F lk (5.10) 

EA l = -iV^djkipj T *? ] )A k (5.11) 

EA l = 4imF 0i + iV^ajkipj ± if j )A k (5.12) 

EF° l = -2imA l + iV2(p, + if l )(f> (5.13) 

EF jk = ±=e ijk ($ i -ir i )4> (5.14) 

Let us mention that for another type of Dirac oscillator- like interaction ~ (ai — 3%) BY '5 the only changes 
are the sign changes at the term if in the fourth and fifth equations of the above system, Eqs. (5.11) 
and (5.12). The two-body Dirac oscillator equations in the form (5.8)-(5.14) could be uncoupled on the 
set containing only functions cj>, <f> and A^ and the another one containing only A^ and F^ v : 

(E 2 - 8m 2 )0 = m - + ir«)<p I I ( 16m ^ f %^ | ^ 

(E 2 - 4m 2 )4> = 2(pi + f l )(p t - if l )4> (5.15) 

EA Q = -2m4> (5.16) 
(E 2 - Sm 2 )A l = 2(pj =f if 3 )(p l ± if v )A 3 - 2(pj =f ir j ){pj ± if^A 1 + 
1 ( /Wmeijkfipf. 

Ml 

and 



+ e\{ (5 - 17) 



EA = (5.18) 
(E 2 - 8m 2 )F° l = A{p, + if l ){pj - if j )F oj - 

Tfi Tfl 

-Ai — (pj ± if 3 )(p l T if l )A 3 + 4i—(pj ± if 3 )(Pj T if J )A l 
E E 

E 2 A l = 2{p 3 ± if 3 )(p l T ir l )A j - 2(pj ± if j )(p 3 T ir^A* + AimEF ' 1 

(5.19) 

(E 2 - Am 2 )F 3k = e ijk e lmn (pi - if l ){pi + if l )F mn . (5.20) 
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This fact proves the Dirac oscillator interaction, like the case of introduction of electrodynamic interaction 
in the Proca or the Bargmann-Wigner equations, does not mix S — 1 and 5 = states. 
Next, the interaction term of the following form: 



v "-?T-flW<*-*> BW (5 ' 21) 

has been deduced |H| from the equation of Relativistic Quantum Constraint Dynamics (RQCR), e.g. 
refs. |52| 33 , or TV- particle Barut equation [Q. In 31 it proved to lead to the Dirac oscillator-like 
interactions provided that the definite choice of the function V(r). In connection with that let us remark 
the curious behavior of the another potential V(r) which has been proposed in ref. [35b, 36]: 

= _ g2 coth(rmn) = _ g2 coth( K r) 
47rr 4nr 

It could be deduced from the one-boson exchange quasipotential V(p, k; E) = —g 2 (p — k)~ 2 by means of 
the transformation into the relatijdstic configurational representation (RCR) using the complete set of 
Shapiro plane-wave functions J37| □: 



(5.24) 



€(j>,r) = (po-prt/m) 1 ~ lrm , p = \/p 2 +m 2 , n = r/\r\. (5.23) 

In the case of the quasipotential ( [5.22| ) the interaction term V, Eq. (5.21), has the different asymptotic 
behavior in three regions (<7 2 /(47r) = 1): 

r{r z — 1) 

(l/r 3 )(ai - a 2 )-Br 5 r, if r >> i and r>l 
-(l/r){d 1 -a 2 )BT 5 r, if i « r < 1, 
behavior in the infrared region (r » -, large distances); and 

V int ~ -2/s(5i - a 2 )Br 5 f, if r « -, (5.25) 

in the ultraviolet region (small distances) . We can convince ourselves that in one of the regions we have 
the Coulomb-like behavior, in the other region, the Dirac oscillator-like behavior. This fact could be 
some quantum field foundations for implementing the Dirac oscillator potential. 

We greatly appreciate valuable discussions with Profs. M. Moshinsky and Yu. F. Smirnov. Authors 
also acknowledge very much Prof. J. S. Dowker for sending reprints of his articles. 



b In the quasipotential approach of Kadyshevsky [ ]35| it is convenient to pass through to the variables Ao = (fcopo — kp) / m 
and A = k(—)p = k~ ^ ^ko — po *jjf m 1 > which have the physical sense of the momentum transfer in the Lobachcvsky space, 
po — p 2 = m 2 . Then, the transformation to the RCR is carried out by means of £(A,r). 
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